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We consider the interactions of a strong gravitational wave with electromagnetic fields using 
the 1+3 orthonormal tetrad formalism. A general system of equations are derived, describing the 
influence of a plane fronted parallel (pp) gravitational wave on a cold relativistic multi - component 
plasma. We focus our attention on phenomena that are induced by terms that are higher order in 
the gravitational wave amplitude. In particular, it is shown that parametric excitations of plasma 
oscillations takes place, due to higher order gravitational non-linearities. The implications of the 
results are discussed. 



I. INTRODUCTION 



There have been numerous investigations on the scattering of electromagnetic waves off gravitational fields (see Refs. 
Previous research has mostly directed its interest towards the effects on vacuum electromagnetic fields (although 
there are exceptions, see e.g. Refs. where the effects of plasmas have been taken into account). Similarly, much 

work concerning gravitational waves have considered the linearized theory, which is obviously the relevant regime for 
gravitational wave detectors, or, in general, for distances far away from the gravitational wave source. Alternatively, 
there has been an interest in exact solutions, and thus a number of exact gravitational wave solutions (see e.g. Ref. 
H and references therein) have been found. In the present paper we will choose an intermediate approach, starting 
with an exact gravitational wave solution, but focusing on a weak amplitude (but still non -linear) approximation, 
and studying the effects induced in a plasma. 

The question under study in this paper is whether non - linear gravitational wave effects - that may be of significance 
close to the gravitational wave source - can give rise to qualitatively new phenomena in plasmas that are absent 
in linearized theory. Close to the source, additional effects apart from non-linearities - due to for example the 
three dimensional geometry and/or the non - radiative part of the gravitational field - are likely to be important 
for astrophysical applications. However, in order to focus on the processes directly induced by non-linearities, a 
somewhat simpler model problem with a unidirectional gravitational wave will be studied: To facilitate the analysis 
of the non -linear interaction between a plasma and a gravitational wave, we make use of the pp-wave solution of 
Einstein's field equations. Furthermore, we introduce a Lorentz tetrad in order to define physical variables in a 
straightforward manner. With this setup, the governing plasma equations can be written in a simple 3 - dimensional 
form. In Maxwell's equations, the gravitational effects are given by effective charge- and current densities. Moreover, 
the fluid equations are given for a cold plasma. 

Previously, parametric excitation of a Langmuir wave and an electromagnetic wave by a linearized gravitational 
wave has been considered ft. Here we address the question whether higher order terms in the gravitational wave 
amplitude can result in new effects, using the above mentioned equations for a cold plasma. In order to demonstrate 
the usefulness of our set of equations, we study the stability properties of a plasma in the presence of a pp- wave. We 
show that including 2 nd order gravitational wave effects may give rise to new phenomena. In particular it is found 
that electrostatic waves can be excited at a resonant surface where the gravitational wave frequency is equal to the 
local plasma frequency. Our results are summarized and discussed in the last section of the paper. 



*E - mail: gert.brodin@physics.umu.se 

^E- mail: mattias.marklund@physics.umu.se 

^E- mail: peter@vishnu.mth.uct.ac.za 



1 



II. PRELIMINARIES 



A. Equations for a general space-time 

We follow the approach presented in ||] for handling gravitational effects in Maxwell's equations. Suppose an 
observer moves with 4- velocity u a (a = 0, ...,3). This observer will measure the electric and magnetic fields |Q| 

E a = F ab u b , B a = \ e abc F bc , (1) 

respectively, where F a b is the EM field tensor. Here e abc is the volume element on hyper - surfaces orthogonal to u a . 

We denote the fluid velocity V a = (7,7V), where 7 = (1 — v 2 )^ 1 / 2 . Let q be the particle charge and n the proper 
number density. Furthermore, we introduce the orthonormal frame (ONF) {e a , a = 0, ■■■■A], each of which is a linear 
combination of the coordinate derivatives d{ = d/dx l , i.e., e a = e a l di. Using the split (|l|) together with j a = qnV a , 
Maxwell's equations V b F ab = j a , V [a F bc] = read 



V • E = p E + p , (2a) 

V • B = p B , (2b) 
E-V X B = -j E -j , (2c) 
B + V X E = -j B , (2d) 

where the "effective" (gravity induced) charge densities and current densities are 

p E = -Y a fja E p - e Q/ ^T%£ 7 , (3a) 

p B = -Y a Pa BP + e a ^T%E 7 , (3b) 

£ = -(rv - r^ )^ + r%E a - e a ^ (r%s 7 + t 6 m b 5 ) , (3c) 

3l = -(r« /9 - v a w)B & + r%i? Q + e at3 ~< (r%i? 7 + r 4 ^) , (3d) 



while p = J2 P s 1^l n anc ^ 3 = J2 P s T1 nv are the matter charge and current densities, respectively (the sums are over 
all particle species). Here Y a bc are the Ricci rotation coefficients with respect to the ONF {e a }. We have introduced 
the notation E = (E a ) — [E X ,E 2 , E 3 ) etc., V = (ei, e2, ea), and the overdot stands for derivative in the direction of 
the time-like vector e . The dot- and cross-products are defined in the usual Euclidean way. 

The energy - momentum tensor for each particle species is assumed to take the form of pressure free matter (dust) , 
T ab = mnV a V b , where to is the rest mass of the particles. Then the conservation equations V b T ab = qnF ab V b give 

e ( 7 n) + V-( 7 w) = - 7 n (rV + T% v a + r%/) , (4a) 
(e + v ■ V) jv = ^{E + v x B) - 7 [r« 00 + (I> + I^o) ^ + T^ 7 ^^] e a . (4b) 



B. Basic relations in the field of a pp - wave 



Previous examinations of interactions between gravitational radiation and EM waves have focused on linearized 
gravitation. On the other hand, one may suspect that there will be interesting effects in the non- linear regime, not 
present to linear order. Below we will show that this is indeed the case. 

In order to address the issue of how strong gravitational radiation may be involved in generation of EM waves, we 
look at the plane fronted parallel (pp) waves (for a discussion, see Refs. ||), in the special case of a linearly polarized 
plane wave 

ds 2 = -dt 2 + a(u) 2 dx 2 + b(u) 2 dy 2 + dz 2 , (5) 

where u — z — t, and a and b satisfy ab uu + a uu b — 0, and the subscript u denotes a derivative with respect to retarded 
time. Note that we have chosen a vacuum geometry, i.e. we have omitted the influence of the plasma on the metric. 
In order to make interpretations simple, we introduce the canonical Lorentz frame 

e = d t , ei = a^d x , e 2 = b^ 1 d v , e 3 = d z . (6) 
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With this frame, the effective charge and current densities (H) read 

p E = -(In ab) u E 3 , (7a) 

p B = -(\nab) u B 3 , (7b) 

j E = -(\nb) u (E 1 -B 2 )e 1 -(lna) u (E 2 +B 1 )e 2 -(\nab) u E 3 e 3 , (7c) 

j B = -{lnb) u (E 2 + B 1 )e 1 + (liia) u (E 1 -B 2 )e 2 -(lnab) u B 3 e 3 . (7d) 

Apart from Maxwell's equations (pa|)-(pd|) [together with the effective charge and current densities (|7a|)-(|7c|)l we 
also need the fluid equations. From the conservation equations (f|) we obtain, the fluid equations using the frame (||): 

d 

— (771) + V-(7ra>) = -fn{\nab) u (l - v\\) , (8a) 
— +vV\jv=^(E + vX B)+7[(lno)„«iei + {lnb) u v 2 e 2 } (1 - v„ ) + 7 [(In a) u u? + (ln&)„uf] e 3 , (8b) 

where U|| = is the velocity parallel to the gravitational wave propagation direction. These equations should be 
satisfied for each particle species. In the limit of small gravitational wave amplitudes and non - relativistic velocities, 
Eqs. together with Maxwell's equations, were given in Ref. ||. All terms with factors (lna&) u are however 

new, and - as we will demonstrate in the remainder of this article - they may induce new phenomena, compared to 
the linear regime. 



III. AN EXAMPLE: PARAMETRIC EXCITATION OF PLASMA OSCILLATIONS 

The longitudinal "currents" and "charges" are second order in the gravitational wave amplitude (see Appendix 
for further details). These second order terms can give rise to qualitatively new phenomena compared to the linear 
regime, and we demonstrate this by considering a simple, but illustrative, example. In what follows, we will investigate 
longitudinal perturbations, i.e., E — (0,0, E), v — (0,0, v) etc., around a cold one-component equilibrium plasma. 
Compared to the case of weak gravitational waves ||, we now have p E B different from zero, and we also have a 
longitudinal contribution to the effective currents. This means that longitudinal EM- and plasma waves can be 
excited. 

In the unperturbed plasma, dn a /dt = 0, E a = 0, and B = |J. We denote the number density perturbation 
by n, i.e., n(z,t) = rip (z) + n(z,t) and assume that all perturbed quantities only depend on t and z. To first order, 
Maxwell's equation ( |2q ) becomes 

dE 

— = (lnab) u E - p qn a v , (9) 



where we have used j m = qnov. Furthermore, the momentum equation (|8b|) becomes 



dv a „ , 

"777 = E . 10 

at m 



Taking the time derivative of Eq. (|9j) and using Eq. (|10|) , we obtain 

^ +u ? p ( z )E=^[{\nab) u E] , (11) 

where u p (z) — [n (z)q 2 po/m] 1 / 2 is the local plasma frequency. Thus the left hand side is the usual equation for 
plasma oscillations in a cold inhomogeneous plasma, and the right hand side is the modification induced by the the 
pp-wave. We next focus ourselves on weak periodic deviations from flat space-time (see the Appendix), where the 
periodicity is 2n/ui. At the resonant surface where u p (z ICS ) — u, we can then have parametric excitation of plasma 
oscillations. We let E(z xes ,t) = E(t) exp [— iut] + c.c, where c.c. denotes the complex conjugate, and assume that 



\dE(t)/dt\ <C uj\E(t)\. At the resonant surface Eq. (|ll| ) then reduces to 



dE 

— = -±iexp(2iiJz rcs )Ljh 2 E* , (12) 
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where the star denotes complex conjugate (see the Appendix). Taking the time-derivative of Eq. ( |l^ ) and using the 
complex conjugate of the same equation, we find E cx exp(ri) where the growth rate is 



Y=\u:\h 2 \. (13) 
Note that the threshold value for excitation is zero, since we have not included any dissipation mechanism of the 



plasma oscillations. Adding electron-ion collision in Eq. (8b), the threshold value /ithr of this instability is of the 
order (^c-j/ojp) 1 / 2 , where v c -\ is the electron-ion collision frequency. 

Clearly, our instability does not occur unless higher order gravitational perturbations are included, in contrast to 
the results in Ref. J(| . Thus the corresponding growth rate is smaller in our case for a given source of gravitational 
radiation. There are still two interesting properties of the above instability as compared to the process in Ref. [[|, 
where parametric excitation of a Langmuir wave and an electromagnetic wave was considered: 

(i) The frequency matching condition in our case is w = uj p , which requires a rather high gravitational frequency 



but is less severe than the condition in Ref. |p||, where w > 2ui p . 



(ii) In contrast to most parametric instabilities in plasmas we have no wave vector matching condition, but instead 
the process takes place at a localized resonance surface z = z Tes where w — uj p (z ics ). This means that there is no 
threshold value for the instability introduced by plasma inhomogeneities. Normally the threshold value is inversely 
proportional to the inhomogeneity scale length fUJI , and close to a binary system, where the effects of gravitational 
radiation are likely to be most important, such a condition for parametric excitation may thus be rather severe. 
Unfortunately, the result of the "no inhomogeneity threshold" depends on the cold plasma approximation, and a finite 
temperature is likely to change the picture. 



IV. SUMMARY AND DISCUSSION 



In the present paper we have investigated a higher order effect of gravitational waves on a plasma. For this purpose 
we have developed a Lorentz tetrad formalism for a cold plasma in the presence of a strong gravitational wave. The 
Lorentz tetrad approach has of course been widely used before, perhaps most notably in the membrane paradigm 
approach to black hole spacetimes [9. The obvious advantage of using a Lorentz tetrad is its direct connection to 
measurements. It is possible to formulate Maxwell's equations such that the gravitational contributions takes the 
form of "charge"- and "current" densities. Similarly, the fluid equations are modified by effective particle sources and 
gravitational forces. Of course, this is not the physical picture behind the equations, but it still provides a useful tool 
for predicting the consequences of the gravitational influence. 

The main purposes of this study has been to (i) provide a framework for investigating strong gravitational pulse 
effects in cold multi - component plasmas, and (ii) to show that higher order gravitational wave effects may be of 
importance, since they introduce effective charges and longitudinal currents, as well as effective "particle sources" 
and gravitational forces. As demonstrated, this in turn makes new processes - such as parametric generation of 
electrostatic waves - possible. Since the effect under discussion is of order h 2 , we do not believe that it will be of 
significance concerning direct earth based observations of gravitational waves. It is possible however that there exists 
favorable circumstances, e.g. close to a binary merger, for which the higher order gravitational effective charge- and 
current densities can play an important role. Close to such sources, the gravitational wave amplitudes can reach 
considerable strength, implying observational possibilities for the induced phenomena. 
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APPENDIX: PERTURB ATIVE EXPANSION OF THE PP - WAVE 



In many situations of interest the gravitational wave amplitude is small, i.e. \a — 1| -C 1, \b — 1| <C 1 and it is 
appropriate to make approximations for the factors (In a) u , (In b) u and (In ab) u that appears in the gravitational source 
terms in Eqs. (0) and (||). We will concentrate on approximately periodic gravitational waves, such as those generated 
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by binary systems, in order to get definite results. Let a(u) = J^^L-oo "™ exp(inwu), b(u) — X)^L-oo ^™ exp(inu>u), 
where a n , b n <C 1, \a n \ ~ |6„| ~ |ai|'"', V n. Furthermore a* = a_„ and similarly for 6. Then a uu b + ab uu — becomes 

oc oo 

(n 2 + TO 2 )a„6 m exp[iw(rt + m)] = . (Al) 



n— — oo m— — oo 



To order, we assume that a = 1 = b Q . To first order, the solution to Eq. (Al) is di = — b\ = h. Clearly, quadratic 



non- linear terms will generate second harmonic terms proportional to exp(2kju). Separating the frequencies in Eq. 



(Al), and concentrating on the second harmonic part we obtain 

2b 2 + 2a 2 - h 2 = . (A2) 

The canonical choice is d 2 = b 2 = (l/4)/i 2 . Physically this means that we minimize the (pseudo) energy density at 
the second harmonic frequency. Thus - for this choice - all the oscillations at 2w are strictly due to the non - linearity 
of Einstein's equations, and no harmonics are assumed to be initially present, i.e. generated by a varying octopole 
moment of the binary source. For astrophysical applications this is not necessarily the most accurate choice (since 
binary systems may indeed have finite octopole moments) , but it has the advantage of clearly isolating the effects due 
to non-linearities. Furthermore, it turns out that including the effect of higher moments of the gravitational source 
(i.e. octupole moments and higher) do not influence our calculations in Sec. III. The reason is that an alternative 



solution to Eq. (A2) ~ a% = (l/4)/i 2 + 5a 2 , b 2 = (l/4)/i 2 — 60,2, instead of b 2 = a 2 = (1/A)h 2 - does not significantly 
affect the factor ln(a6) u , since it is independent of Sa 2 to second order in the gravitational amplitude, provided 
Sa 2 ~ h 2 . 

Continuing to 3 rd order, a similar calculation shows that we can make the natural choice 0,3 = 63 = 0. For the 4 th 



order terms, Eq. ( |Al| ) gives 

32(d 4 + 64) - h 4 = . (A3) 

where the canonical choice d 4 = 64 = (l/64)/i 4 is made. Continuing this procedure, it turns out that all terms odd in 
n 7^ 1 disappear, while the terms even in n satisfy a n = b n V n. 

Using the above results, the logarithmic factors in Eqs. (0) and @ become 

(lna) u = iuihexp(iuju) — ^iojh 2 exp(2ia;tt) + jiujh 3 exp(3iww) — j^^h 4 exp(4ia;ii) + c.c (A4a) 
(lnb) u = —iojhexp(iuu) — ^iuih 2 exp(2ia->w) — jiuih 3 exp(3kju) — j^^h 41 exp(4iww) + c.c (A4b) 
(lna6) u = —iwh 2 exp(2iuju) — exp(4iam) + c.c (A4c) 

to 4 th order in the gravitational amplitude. This procedure may of course be continued to arbitrary order, noting 
that this in general will result in an asymptotic series, i.e., it does not necessarily converge towards a solution of 

a U ub + ab uu = 0. 
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